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Axial  waves  predicted  by  a  two-dimensional  hybrid  numerical  model  have  been  used  to 
estimate  electron  cross  field  transport  due  to  tilted  waves  with  azimuthal  components.  Since 
the  radial-axial  hybrid  simulation  cannot  model  these  tilted  waves  directly,  the  predicted 
axial  waves  are  assumed  to  couple  symmetrically  into  two  counter-propagating  axial- 
azimuthal  waves.  A  linearized  two-dimensional  dispersion  relation  is  solved  to  obtain  the 
azimuthal  component  of  wavenumber  consistent  with  the  frequency  and  axial  wavenumbers 
predicted  by  the  radial-axial  simulation.  The  resulting  transport  profiles  are  in  qualitative 
agreement  with  experimental  measurements  of  electron  mobility  when  the  power  levels  of 
the  fluctuations  approach  saturation  levels.  For  large  amplitude  fluctuations,  the  predicted 
electron  mobility  profile  exhibits  a  transport  barrier  near  the  exit  plane  of  the  thruster 
consistent  with  experimental  measurements,  suggesting  that  this  barrier  may  be  due  to  a 
transition  from  transport-enhancing  azimuthal  waves  to  transport-limiting  axial  waves  in 
this  region. 
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Nomenclature 

magnetic  field,  T 

electric  field,  V/m 

current  density.  A/m2 

temperature,  K 

velocity,  m/s 

electron  charge,  C 

arbitrary  fluctuating  quantity 

Boltzmann  constant,  J/K 
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k 

= 

perturbation  wavenumber,  1/m 

m 

= 

mass,  kg 

n 

= 

number  density,  1/m3 

P 

= 

pressure,  N/m2 

<7 

= 

charge,  C 

t 

= 

time,  s 

X 

= 

position  vector 

r 

= 

radial  direction 

z 

= 

axial  direction 

e 

= 

azimuthal  direction 

a 

= 

ionization  rate,  1/s 

P 

= 

mobility,  m2/Vs 

V 

= 

collision  frequency,  1/s 

<t> 

= 

electric  potential,  V 

CO 

= 

perturbation  frequency,  rad/s 

C0c 

= 

cyclotron  frequency,  rad/s 

(OJ 

= 

Hall  parameter 

Subscript 

e 

= 

electron 

i 

= 

ion 

o 

= 

time  average 

r 

= 

radial  direction 

z 

= 

axial  direction 

e 

= 

azimuthal  direction 

I.  Introduction 

A  Hall  thruster  is  an  electric  propulsion  device  capable  of  providing  continuous  low  thrust  while  maintaining 
high  efficiency  and  propellant  utilization. 1  In  a  typical  Hall  thruster,  electrons  travel  from  an  external  cathode 
to  an  anode  at  the  base  of  an  annular  channel.  Along  their  path,  the  electrons  encounter  a  region  of  high  resistivity 
due  to  an  imposed  radial  magnetic  field,  and  become  trapped  in  an  azimuthal  E  x  B  Hall  current.  When  neutral 
particles  injected  from  the  anode  encounter  this  region  of  high  electron  density,  electron  impact  ionization  occurs. 
The  newly  born  unmagnetized  ions  then  accelerate  out  of  the  thruster  due  to  the  electrostatic  potential  maintained 
between  the  anode  and  cathode.  While  this  process  results  in  a  high  specific  impulse  device  suitable  for  various  low 
thrust  applications,  the  physics  behind  Hall  thruster  operation  is  not  well  understood,  particularly  the  mechanism 
which  allows  electrons  to  cross  magnetic  field  lines. 

Electron  cross-field  transport  in  Hall  thrusters  is  experimentally  observed  to  be  higher  than  predicted  by  classical 
theory.2  One  possible  mechanism  producing  anomalously  high  electron  current  is  correlated  azimuthal  fluctuations 
in  the  electric  field  and  plasma  density  within  these  devices.  The  most  widely-used  numerical  model  for  predicting 
Hall  thruster  behavior  is  a  radial-axial  hybrid  fluid-particle  code.  While  this  model  is  capable  of  producing  many  of 
the  features  observed  within  Hall-effect  thrusters,  such  as  the  commonly  observed  "breathing  mode"  oscillations,3 
models  of  this  type  require  empirical  fitting  parameters  to  accurately  predict  electron  cross-field  mobility.  Since 
radial-axial  hybrid  models  assume  axisymmetric  plasma  properties,  they  are  incapable  of  capturing  the  enhanced 
transport  that  results  from  azimuthal  fluctuations.  However,  while  these  models  lack  resolution  in  the  Hall  direction, 
they  generate  a  wide  spectrum  of  frequencies  and  wavelengths  corresponding  to  axial  waves.4  Although  the  relative 
phase  of  density  and  velocity  fluctuations  in  the  axial  direction  does  not  produce  anomalous  transport,  these  axial 
waves  may  couple  into  tilted  axial-azimuthal  waves  which  are  capable  of  explaining  the  observed  electron  mobility. 

This  work  will  use  a  small-amplitude  linear  perturbation  analysis  to  model  the  coupling  of  the  axial  waves 
generated  by  the  radial-axial  hybrid  simulation  into  tilted  axial-azimuthal  waves.  Despite  the  distribution  of 
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simulated  power  in  the  hybrid  Hall  thruster  model  over  a  variety  of  wavelengths  and  frequencies,  the  analysis  will 
be  simplified  by  considering  only  the  dominant  frequency  at  a  variety  of  wavelengths.  The  computed  spatially- 
varying  electron  cross  field  transport  within  the  interior  of  the  Hall  thruster  channel  caused  by  the  inferred  axial- 
azimuthal  waves  will  be  presented.  Despite  using  a  linear  analysis  to  predict  behavior  in  the  non-linear  saturation 
regime,  the  computational  results  will  demonstrate  that  correlated  fluctuations  are  capable  of  enhancing  electron 
transport  in  Hall  thrusters.  Additionally,  results  of  the  linear  analysis  suggest  that  the  transport  barrier 
experimentally  observed  near  the  thruster  exit  may  be  due  to  a  transition  from  transport -enhancing  azimuthal  waves 
to  non-enhancing  axial  waves. 


II.  Numerical  Models 


A.  Hybrid  Hall  Thruster  Model 

The  numerical  simulation  used  to  calculate  axial  wave  properties  is  a  radial-axial  hybrid  fluid  particle-in-cell 
(PIC)  model5  6  based  on  work  by  Fife.7  The  model  employs  a  quasi-one-dimensional  fluid  treatment  of  electrons  and 
a  three-dimensional  PIC  treatment  of  the  heavy  species,  Xe  and  Xe+.  The  two  solutions  are  coupled  assuming  space- 
charge  neutrality. 

The  computational  geometry  used  in  the  simulation  corresponds  to  the  laboratory  discharge  referred  to  here  as 
the  Stanford  Hall  thruster.s  The  annular  channel  is  approximately  8  cm  in  length  and  1.2  cm  in  width.  The 
computational  grid,  which  includes  a  two-dimensional  slice  of  the  channel  and  near  field  region,  is  shown  in  Fig.  1. 
The  insulating  channel  walls  are  made  of  alumina.  A  mass  flow  rate  of  2  mg/s  is  implemented  in  the  model  to  match 
experimental  conditions.  Measurements  of  the  magnetic  strength  in  the  radial  and  axial  directions  along  the  channel 
centerline  are  used  to  impose  a  constant  external  magnetic  field.  A  discharge  voltage  of  200  V  is  considered  in  the 
results  presented  here. 

The  electron  fluid  is  governed  by  a  continuity  equation,  momentum  equations  parallel  and  perpendicular  to 
magnetic  field  lines,  a  one -dimensional  electron  energy  equation,  and  a  current  continuity  equation.  Perpendicular  to 
magnetic  contours,  an  experimentally-based  electron  cross  field  mobility9  is  imposed  in  order  to  calculate  the 
electron  cross  field  velocity  and  thermal  diffusivity.  The  electron  mobility  along  magnetic  field  lines  is  assumed  to 
be  infinite  resulting  in  isothermal  magnetic  field  contours.  Therefore,  along  field  lines,  the  electric  force  is  balanced 
by  a  pressure  force  resulting  in  a  Boltzmann  description  of  the  electrons.  The  electron  energy  equation  includes  a 
Joule  heating  energy  source  term  and  ionization  and  wall  damping  energy  sink  terms.  The  wall  damping  treatment  is 
similar  to  that  of  Barral  et «/. 10  However,  while  Barral  calculates  a  separate  electron  temperature  parallel  to  magnetic 
contours  in  order  to  determine  the  heat  transfer  to  the  wall,  this  model  lowers  the  effective  temperature  at  the  wall  to 
account  for  anisotropy  and  the  non-Maxwellian  nature  of  the  electron  velocity  distribution  function.  The  effective 
temperature  used  to  compute  wall  energy  loss  producing  the  best  agreement  with  experiment  was  found  to  be  40% 
of  the  perpendicular  temperature  computed  from 
the  ID  energy  equation. 

The  motion  of  the  heavy  species,  Xe  and  Xe+,  is 
solved  in  three  dimensions  using  cylindrical 
coordinates.  Neutral  particles  are  injected  at  the 
anode  and  scatter  off  channel  walls  assuming  a  one¬ 
way  Maxwellian  flux  distribution.  Ions  which 
impact  channel  walls  are  neutralized  before  being 
re-emitted  into  the  channel.  Ionization11'14  and 
charge  exchange  collisions15  take  place  everywhere 
inside  the  computational  domain  based  on  local 
plasma  properties.  Neutral  particles  are  injected 
from  the  computational  boundaries  in  the  near-field 
region  in  order  to  simulate  a  back  pressure  of  0.05 
mTorr. 

For  computational  manageability, 

macroparticles  are  used  to  represent  large  groups  of 
neutrals  and  ions  rather  than  simulating  individual 
particles.  Since  neutral  and  ion  densities  differ  by 
orders  of  magnitude  over  the  length  of  the  domain.  Figure  1.  Computational  grid  of  the  Stanford  Hall 

thruster  radial-axial  simulation. 
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the  size  of  the  superparticles  vary  with  both  space  and  species.  The  simulation  is  initialized  with  approximately 
500,000  marcroparticles  of  each  heavy  particle  species.  The  ion  and  neutral  timestep  is  25  ns,  while  the  electron 
timestep  is  0.1  ns.  The  simulation  completes  625  ps  in  approximately  one  day  on  a  3.8  GHz  Pentium  4  processor. 

B.  Linearized  Axial-Azimuthal  Wave  Model 

As  shown  by  Fernandez  et  al.,  the  axial  waves  generated  by  a  radial-axial  hybrid  Hall  thruster  simulation  can  be 
understood  using  a  simple  linearized  perturbation  analysis.4  However,  due  to  the  relative  phase  of  the  density  and 
velocity  fluctuations,  the  axial  waves  predicted  by  the  linear  analysis  cannot  be  responsible  for  the  observed 
anomalous  electron  cross  field  transport.  Also,  laboratory  measurements  may  suggest  that  these  waves  do  not  exist 
in  reality.  While  it  is  possible  that  the  physics  included  in  the  radial-axial  simulation  are  insufficient  to  describe 
experimental  trends,  another  explanation  is  that  the  power  produced  in  these  axial  waves  is  dissipated  through  wave 
coupling.  Experimentally,  physical  grooves  in  the  ceramic  channel  are  found  to  develop  at  the  exit  plane  of 
laboratory  Hall  thrusters  after  many  hours  of  operation,  suggesting  the  presence  of  standing  waves.  Therefore,  in  the 
analysis  presented  here,  it  is  assumed  that  the  axial  waves  generated  by  the  radial-axial  simulation  couple 
symmetrically  into  two  counter-propagating  tilted  waves  with  half  the  frequency  and  half  the  axial  wavenumber,  as 
shown  in  Fig.  2.  Using  a  similar  linear  analysis  as  described  in  Ref.  4,  a  two-dimensional  linearized  dispersion 
relation  can  be  derived  in  order  to  compute  the  azimuthal  components  of  wavenumber  consistent  with  the  driving 
axial  wave. 

In  this  analysis,  a  fluid  model  is  used  to  describe  the  electrons  and  ions,  and  a  local  assumption  is  enforced 
which  neglects  spatial  non-uniformity  of  steady  state  properties.  Small  amplitude  perturbations  are  assumed  in  order 
to  linearize  the  continuity  and  momentum  equations  for  the  electrons  and  ions.  All  densities,  n,  and  velocities,  U,  as 
well  as  the  electric  potential,  <j>,  are  assumed  to  be  composed  of  a  steady  state  component,  f  ,  and  a  fluctuating 

component,  /  exp(  i(k  ■  x  -  cot)) ,  as  shown  below. 

/  =  f„  +  f  expW^  •  v  -  ox))  (1) 

In  this  equation,  k  is  the  wavevector  and  co  is  the  frequency  of  the  perturbation.  Note  that  the  fluctuation  amplitude, 
f  ,  may  be  complex  in  order  to  account  for  differences  in  phase,  x  represents  an  arbitrary  position  vector  such  that: 

k  -x  =  k.z  +  kgd  (2) 

Despite  the  cylindrical  nature  of  the  problem,  the  coordinates  z  ,  0 ,  and  >■  are  treated  as  Cartesian.  Note  that  &  is 
assumed  to  be  zero. 

In  the  electron  momentum  equations,  inertia  terms  are  retained  and  a  classical  description  based  on  electron- 
neutral  collisions  is  used  to  describe  the  electron  cross  field  mobility.  Due  to  the  predominantly  axial  steady  state 
electric  field,  the  steady  state  ion  velocity  is  assumed  to  be  purely  axial.  However,  the  small-amplitude  fluctuating 
component  of  ion  velocity  will  retain  contributions  in  both  the  axial  and  azimuthal  directions.  Rather  than 
incorporating  a  linearized  version  of  the  electron  energy  equation,  the  electron  temperature  is  assumed  not  to  exhibit 
oscillatory  behavior.  Lastly,  the  system  of  equations  is  closed  assuming  space-charge  neutrality. 

The  system  of  equations  to  be  solved  is  as  follows: 

1.  Electron  Continuity  Equation 

c\  +  aMJ+4c/j=m  (3) 

dt  dz  dd 


2.  Ion  Continuity  Equation 


dn,  |  d(n,Uiz)  |  d(n,U  l0) 
dt  dz  dO 


=  anj 


(4) 


3. 


4. 


Electron  Axial  Momentum  Equation 

dU  dU  dU 

+  £/  — ^  +  U, 


dt  ~ e!:  dz  ~ de 

Electron  Azimuthal  Momentum  Equation 
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Figure  2.  Illustration  of  axial  driving  wave  with  given  axial  wavenumber  and  frequency  coupling  into  two 
symmetric  axial-azimuthal  waves  with  half  the  axial  wavenumber  and  frequency. 


5.  Ion  Axial  Momentum  Equation 


dU  dU .  e 

” -  +  U  „ — ^  =  — -aUh 


dt 


dz  m. 


(7) 


6.  Ion  Azimuthal  Momentum  Equation 


dU.„  8U.„  e 

Ee-aUie 


dt 


dz  m. 


(8) 


7. 

Equation  of  Quasineutrality 

ne  =  «, 

(9) 

8. 

Equation  for  an  Electrostatic  Field 

E  =  -V  <j> 

(10) 

9. 

Equation  for  Electron  Pressure 

p  =  n  k.T 

re  e  b  e 

(11) 

A  detailed  description  of  the  solution  of  this  system  of  equations  through  substitution  of  quantities  of  the  form 
given  by  Eq.  (1)  is  provided  in  Chapter  8  of  Ref.  6.  The  linear  analysis  provides  equations  for  the  fluctuation 
amplitudes  in  terms  of  steady  state  quantities  and  wave  properties.  A  summary  of  the  resulting  equations  is  given 
below: 

1.  Electron  Density  Fluctuation  Amplitude 


2. 


3. 


m 

(co  -  k  ■  Ueo  +  i  ve )~  -  co2e 

/■  _  _ 

\ 

co-k  -  Uei 

,  -  ia 

k 2 

^co-k-  Ueo 

+  ive  j 

me 

Ion  Density  Fluctuation  Amplitude 


k2</ 


i  Uo  I  ,  Tt  \2  2 

m,  \a)-kzUko)  +a- 

Axial  Electron  Velocity  Fluctuation  Amplitude 


U._  = 


-  e 

- 1 

m„ 


kje 


[co-k-Ueo+ive\ -ml 


-((a-k-Ueo+ive)k:  +  icocekg ) 


(12) 


(13) 


(14) 
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4.  Azimuthal  Electron  Velocity  Fluctuation  Amplitude 
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5.  Axial  Ion  Velocity  Fluctuation  Amplitude 
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6.  Azimuthal  Ion  Velocity  Fluctuation  Amplitude 
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In  the  above  equations,  the  cyclotron  frequency  is  defined  as: 

c/B 

co„  = 


(15) 


(16) 


(17) 


(18) 


Once  the  fluctuation  amplitudes  have  been  found,  a  dispersion  relation  can  be  derived  using  the  quasineutral 
assumption  (Eq.  (9)).  Setting  the  electron  density  perturbation  amplitude  equal  to  the  ion  density  perturbation 
amplitude  yields: 


e{[co-k-Ueo+ive) 
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This  equation  can  be  further  simplified  by  defining  the  dot  products  as: 

k-Ue.0=kUMp 
k  ■  Ui  o  =  kzUiz  o  =  kUalp 


(19) 

(20) 
(21) 

where  the  components  of  velocity  parallel  to  the  wave  traveling  at  angle  y  with  respect  to  the  ’-direction  are  given 
by: 

=  Uezo  cos  y  +  Ue6o  sin  y  (22) 

Um„=Uizocosy  (23) 

Using  these  substitutions  and  expanding  Eq.  (18)  in  terms  of  the  angular  frequency  yields  the  following 
simplified  dispersion  relation: 


where: 


C3a)  +  C2co  +C,j{y  +  C0  —0 


C,  =  m,  +  m. 


(24) 

(25) 

(26) 
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C2  =-(2me+ml\kUMio-ive)-me(kUMp+ia)-2inikUmp 
C,  =  me((kUelu>  -ivej  -  co2e)~  khTek2  +m,(k2U2p -a2).. 

....+  2mi(kUMp  —  ive)kUmp  +  2me(kUeHp -ive\kUeHp +ia) 

C0  =mc-(kUenp-ive)2(kUellp+ia)+me(kUenp+ia}.. 

■■■■-mi(kUeUp -ive)k2U2p  +  (kUMfi - ive)kbT,k 1  +m{kUMe - ive)a 2 

Given  the  frequency  of  the  perturbation  and  the  axial  component  of  wavenumber,  azimuthal  components  of 
wavenumber  can  be  found  which  are  consistent  with  this  dispersion  relation. 

C.  Linearized  Fluctuation-Induced  Mobility  Model 

In  the  radial-axial  hybrid  simulation,  since  azimuthal  forces  are  neglected  due  to  axisymmetry,  the  axial  electron 
current  density  is  related  to  the  axial  forces  through  the  electron  cross  field  mobility,  u .  Assuming  that  the  magnetic 
field  is  primarily  radial  such  that  the  perpendicular  direction  can  be  expressed  by,  z,  the  axial  electron  current 
density  is  given  by: 


Jez  =  eneii\ 


E_  + 


1  dpe 
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Classically,  the  mobility  is  given  by: 


(30) 


In  most  Hall  thrusters,  the  electron  cyclotron  frequency,  coce,  dominates  the  electron  collision  frequency  with 
heavy  particles,  y  ,  in  all  regions  of  the  channel  such  that  the  cross  field  mobility  can  be  expressed  as: 


M  = 


1 


(31) 


where  ®  z  is  the  dimensionless  quantity  called  the  Hall  parameter.  Neglecting  azimuthal  forces,  the  azimuthal 
electron  current  density  is  given  by: 

1  dp. 
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Therefore,  an  effective  Hall  parameter  can  be  defined  in  terms  of  the  azimuthal  and  axial  electron  current  densities 
as  follows: 


(c 


r) 

‘’eg  , 

J  0-7 


(33) 


In  the  radial  axial  hybrid  simulation  used  to  generate  the  axial  waves  on  which  this  analysis  is  based,  the 
effective  Hall  parameter  was  measured  experimentally.2  In  many  numerical  models,  the  effective  mobility  is  treated 
using  adjustable  fitting  parameters  which  are  optimized  to  produce  agreement  with  experimental  measurements. 
However,  in  order  to  simulate  new  thrusters  or  existing  thrusters  at  operating  conditions  other  than  those  for  which 
experimental  measurements  have  been  made,  the  mechanism  producing  the  anomalous  transport  should  be 
understood.  In  this  linear  analysis,  an  effective  Hall  parameter  is  calculated  using  simulated  plasma  properties  and 
inferred  axial-azimuthal  wave  properties  which  are  used  to  compute  an  anomalous  axial  electron  current.  The 
fluctuation-based  Hall  parameter  is  therefore  determined  using: 

r  \  ed ) classical  (34) 


co... 
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In  reality,  the  azimuthal  current  density  will  also  consist  of  classical  and  fluctuating  components.  However,  since 
the  radial-axial  hybrid  simulation  assumes  uniformity  in  the  azimuthal  direction,  the  azimuthal  current  density  is 
assumed  in  the  numerical  model  to  be  given  by  Eq.  (32)  where  all  quantities  are  calculated  based  on  time-averaged 
quantities.  Also,  in  the  above  expression. 


the  classical  axial  current  density  is  given 
by  Eq.  (29)  where  in  addition  to  using 
time-averaged  quantities,  the  cross  field 
mobility  is  based  only  on  classical 
collisions  between  electrons  and  heavy 
particles. 

III.  Results  and  Discussion 

A.  Simulated  Axial  Waves 

By  imposing  the  experimentally 
measured  electron  cross-field  mobility 
profile,  the  radial-axial  hybrid  Hall 
thruster  simulation  exhibits  various 
fluctuations  in  addition  to  the 
experimentally  observed  "breathing 
mode"  oscillations.  As  an  illustration.  Fig. 
3  demonstrates  a  simulated  dispersion  plot 
corresponding  to  a  location  near  the  exit 
plane  of  the  thruster  computed  using  a 


Axial  Position:  8  cm 


-400  -200  0  200 

Inverse  Wavelength  [1/m] 


400 


Figure  3.  Dispersion  plot  at  an  axial  location  of  8  cm  along  channel 
centerline,  corresponding  to  the  exit  plane.  The  colorbar  represents 
loglO  of  total  power  in  [V2]. 
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Dominant  Frequency 


Figure  4.  Spatial  dependence  of  simulated  dominant  frequency  for 
three  different  axial  wavenumbers. 


wavelet  analysis  of  the  simulated,  time- 
dependent,  spatially-varying  plasma 
potential.  Due  to  both  the  transient  and 
spatially  varying  nature  of  the  plasma 
properties,  a  Morlet  wavelet  analysis16 
was  selected  as  a  means  to  capture  both 
the  time-dependent  frequency  response, 
as  well  as  reduce  both  noise  and  edge 
effects. 

As  shown  Fig.  3,  the  simulated 
power  is  spread  over  a  wide-range  of 
frequencies  and  wavenumbers.  The 
dominant  frequency  is  found  to 
increase  as  a  function  of  increasing 
wavenumber,  or  decreasing 
wavelength.  As  previously  mentioned, 
Fernandez  et  til.'  has  shown  that  the 
relationship  between  wavenumber  and 
frequency  exhibited  by  the  hybrid 
model  can  be  replicated  through  a 
linear  analysis  of  the  equations  input 


into  the  numerical  model.  In  order  to 
calculate  the  resulting  transport  due  to 
wave-coupling,  a  more  thorough 
analysis  would  account  for  the  wide 
range  of  excited  frequencies  and 
wavenumbers.  However,  in  this 
analysis,  only  a  dominant  frequency 
and  wavenumber  are  considered.  In 
order  to  evaluate  the  influence  of 
different  wavelength  oscillations,  three 
representative  wavenumbers  have  been 
chosen  for  examination:  100,  300,  and 
500  m"1.  Figure  4  shows  the  dominant 
frequency  as  a  function  of  axial 
location  along  the  centerline  for  each 
of  these  three  axial  wavenumbers.  The 
dominant  frequency  was  found  to  vary 
between  ~  500  kHz  and  8  MHz. 

B.  Predicted  Axial-Azimuthal  Wave 
Coupling 

Using  the  axial  wave  properties 
predicted  by  the  radial-axial  hybrid 
simulation,  as  shown  in  Fig.  4,  the 
properties  of  the  corresponding  axial- 
azimuthal  waves  can  be  computed.  As 
discussed  in  Section  IIB,  the  axial 
driving  wave  is  assumed  to  couple 
symmetrically  into  two  counter- 
propagating  axial-azimuthal  waves. 
Therefore,  each  generated  wave  is 
assumed  to  have  half  the  frequency  and 
half  the  axial  wavenumber  as  the 


100 


Axial  Position  [m] 


Figure  5.  Illustration  of  axial  driving  wave  with  given  axial 
wavenumber  and  frequency  coupling  into  two  symmetric  axial- 
azimuthal  waves  with  half  the  axial  wavenumber  and  frequency. 
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simulated  axial  driving  wave.  Using  these  frequencies  and  wavenumbers  at  each  location,  along  with  time-averaged 
properties  computed  by  the  simulation,  the  fluctuation  amplitudes  and  the  azimuthal  components  of  the 
wavenumbers  can  be  computed  through  solution  of  the  dispersion  relation  as  given  by  Eqs.  (24)-(28). 

Figure  5  illustrates  the  computed  magnitude  and  angle  of  propagation  for  the  tilted  waves  for  each  of  the  three 
driving  wave  axial  wavenumbers.  Note  that  more  than  one  magnitude  and  angle  is  possible  at  each  location 
indicating  that  more  than  one  solution  exists  to  the  dispersion  relation  which  is  consistent  with  the  prescribed 
frequency  and  axial  wavenumber.  As  shown  in  the  upper  plot  of  Fig.  5,  the  inverse  wavelength  of  the  tilted  wave  is 
found  to  vary  between  50  and  700  m’1  which  corresponds  to  wavelengths  between  1.5  mm  and  2  cm.  The  lower  plot 
of  Fig.  5  shows  that  in  addition  to  the  predominantly  axial  root  near  zero  degrees,  there  are  two  counter-propagating 
solutions,  as  expected,  with  propagation  angles  approaching  +/-  90  degrees  near  the  anode  corresponding  to 
azimuthal  waves  which  become  more  axial  toward  the  exit  plane  at  8  cm.  The  smaller  wavenumber  case, 
corresponding  to  the  largest  considered  driving  wavelength  of  1  cm,  remains  more  strongly  azimuthal  for  a  greater 
region  of  the  channel. 

C.  Computed  Mobility 

Using  the  wave  properties  predicted  in  the  previous  section,  combined  with  the  steady  state  values  calculated  by 
the  radial-axial  hybrid  model,  an  effective  electron  cross  field  Hall  parameter  can  be  computed  from  Eq.  (34).  In 
addition  to  the  classical  contribution  to  axial  current  based  on  time-averaged  values,  the  computed  Hall  parameter  is 
also  based  on  correlated  fluctuations  between  the  electron  density,  given  by  Eq.  (12),  and  the  axial  electron  velocity, 
given  by  Eq.  (14).  The  one  remaining  unknown  parameter  necessary  in  order  to  calculate  a  mobility  profile  is  the 
power  in  the  fluctuation  related  to  the  fluctuation  amplitude  of  the  potential,  <j>  .  The  fluctuation  amplitude  predicted 
by  the  wavelet  analysis,  represented  by  the  solid  black  curve  in  Fig.  5,  is  insufficient  to  explain  the  anomalous 
electron  transport  observed  in  a  Hall  thruster.  However,  the  wavelet-based  potential  fluctuation  amplitude  is 
calculated  based  on  one  narrow  band  in  frequency  and  wavenumber.  As  illustrated  by  Fig.  3,  power  is  distributed 
over  a  wide  spectrum  of  frequencies  and  wavenumbers  which  have  been  neglected  in  this  analysis.  Therefore,  in 
order  to  estimate  the  possible  effect  of  including  too  little  power,  a  potential  fluctuation  amplitude  has  been 
computed  using  a  saturation  model.  This  model,  which  defines  the  upper  limit  of  the  power  contained  in  the  wave, 
computes  the  potential  fluctuation  amplitude  necessary  such  that  the  electron  density  perturbation  amplitude,  he,  is 

equal  to  the  steady  state  value,  neo,  in  Eq.  (12).  For  comparison,  the  total  power  in  three  different  frequency  bands 


Figure  6.  Comparison  of  potential  fluctuation  amplitude  computed  using  saturation  model 
and  wavelet  analysis  with  total  fluctuation  amplitude  in  various  frequency  bands. 
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- Classical 
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Figure  7.  Computed  inverse  Hall  parameter  using  linear  perturbation  model  assuming  axial  - 
azimuthal  waves  coupled  to  simulated  axial  waves.  Saturation  potential  fluctuation  is  used  rather 
than  simulated  value. 

computed  from  the  simulation  results  is  shown  as  well.  Note  that  in  the  upstream  region  of  the  channel,  the  total 
simulated  power  is  more  than  enough  to  saturate  the  wave.  Downstream,  the  simulated  variation  in  potential  is  not 
large  enough  to  result  in  perturbation  amplitudes  on  the  order  of  the  steady  state  electron  density. 

As  an  upper  limit,  the  saturated  potential  amplitude  has  been  used  to  compute  electron  transport  from  Eq.  (34). 
The  computed  inverse  Hall  parameter  is  compared  with  the  experimentally  measured  inverse  Hall  parameter  in  Fig. 
7.  The  results  indicate  that  one  root  of  the  solution  corresponds  to  classical  transport,  while  a  second  root  results  in 
anomalous  transport  upstream  near  the  anode.  The  region  of  anomalous  transport  extends  farther  into  the  channel  for 
smaller  wavenumbers  (larger  wavelengths).  While  the  model  does  predict  a  small  region  near  7  cm  (approximately 
1  cm  upstream  of  the  exit  plane)  where  the  transport  is  also  anomalously  high,  in  general,  transport  is  seen  to  reduce 
near  the  exit  place  in  the  region  where  a  transport  barrier  is  experimentally  found  to  exist.  In  the  linearized  model, 
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the  region  of  reduced  mobility  is  found  to  correspond  to  a  transition  from  azimuthal  to  axial  waves.  Note  that  the 
best  agreement  occurs  for  the  largest  considered  driving  wavelength  of  1  cm  (100  m'1).  Although  fluctuation 
magnitudes  increase  near  the  exit  plane,  the  relative  phase  of  the  density  and  velocity  perturbations  result  in  less 
overall  transport,  reproducing  experimentally  observed  trends. 

IV.  Summary 

A  linearized  perturbation  analysis  has  been  performed  to  examine  the  possibility  of  wave  coupling  between  the 
axial  waves  generated  by  a  radial-axial  hybrid  Hall  thruster  model  and  axial-azimuthal  waves  which  may  be 
responsible  for  observed  anomalous  transport.  The  predicted  tilted  waves  have  frequencies  on  the  order  of  several 
hundred  kilohertz  near  the  anode  to  a  few  megahertz  near  the  exit  plane.  The  wavelengths  considered  are  on  the 
order  of  millimeters  to  one  centimeter.  The  angles  of  propagation  of  the  tilted  waves  consistent  with  a  dispersion 
relation  derived  based  on  a  linearized  analysis  suggest  that  the  tilted  waves  transition  from  almost  purely  azimuthal 
near  the  anode  to  nearly  axial  at  the  exit  plane.  Using  saturation  amplitudes  for  the  potential  fluctuations,  the 
computed  transport  was  found  to  qualitatively  agree  with  experimental  measurements.  The  predicted  mobility  is 
found  to  be  anomalously  high  near  anode  and  in  the  first  half  of  the  channel  and  transition  to  classical  values  in  the 
second  half  of  the  channel  near  the  exit  plane.  Based  on  the  linear  perturbation  model,  the  reason  for  this  transition 
is  found  to  correspond  to  a  transition  from  transport-enhancing  azimuthal  waves  near  the  anode,  where  the  density 
and  velocity  perturbations  are  in  phase,  to  transport-suppressing  axial  waves  near  the  exit  plane,  where  the 
fluctuations  are  nearly  90  degrees  out  of  phase. 
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